Elasticity-driven interaction between vortices in type-II superconductors 
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The contribution to the vortex lattice energy which is due to the vortex-induced strains is cal- 
culated covering all the magnetic field range which defines the vortex state. This contribution is 
compared with previously reported ones what shows that, in the most part of the vortex state, it 
has been notably underestimated until now. The reason of such underestimation is the assumption 
that only the vortex cores induce strains. In contrast to what is generally assumed, both core and 
non-core regions are important sources of strains in high-K superconductors. 
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I. INTRODUCTION 

Since long time ago, much attention has been paid 
to the role of long-range strain fields in the vor- 
tex state of type-II superconductors. It is well 
known, for instance, that interaction between defect- 
induced strains and vortices causes pinning phenom- 
ena. These phenomena have been extensively stud- 
ied almost since Abrikosov predicted the superconduct- 
ing vortices (see, e.g., Refsi2i2i2i£i£iL£). It is also 
known that vortex-induced strains give a contribution 
to the energies of the vortex lattices (VL). This con- 
tribution proves to be essential when discussing the ob- 
served correlations^ between VL's and crystal lattices in 
anisotropic superconductors The vortex-induced 

strains might be important in vortex inertia also, because 
they contribute to the effective masses of vortices^ 

In this paper, we calculate the contribution to the VL 
energy which is due to the vortex-induced strains. Com- 
parison with the previously reported calculations i2*iii£ 
shows that, for magnetic fields not so close to the up- 
per critical field H C 2, this contribution has been notably 
underestimated until now. The reason of such underesti- 
mation is connected with the fact that, contrary to what 
is assumed in many occasions, the vortex core is not the 
primary source of strain when the Ginzburg-Landau pa- 
rameter k of the superconductor is large. 

To clarify this point we shall revise, first of all, the 
strain induced by a single vortex. This strain is due to all 
the spatial variations of the density of superconducting 
electrons that the vortex provokes. The vortex core is 
a region of strong variations, but is not the only one. 
There also exist a region of smooth variations which is 
associated with the presence of superconducting currents. 
In high-K superconductors, the size of the latter region is 
much larger than that of the core. Just because of this 
greater extension, the non-core variations of the density 
of superconducting electrons are what finally emerge as 
the main sources of strains. 

Previous calculationsifliiiiiS of the elasticity-driven in- 
teraction between vortices was based on models that as- 



sume, from the beginning, that only the vortex cores in- 
duce strains. So an important source of strains in high-K 
superconductors is overlooked in all these works. But 
note that, even doing so, it was shown that this inter- 
action was strong enough to explain the observed corre- 
lations between VL's and crystal lattices in NbSe2. We 
revise this elasticity-driven interaction showing that the 
proper inclusion of all the sources of strain increases its 
importance in the corresponding problems. 

Let us mention that we evaluate this interaction tak- 
ing into account all the elastic degrees of freedom of free 
samples of finite size, i.e. taking into account that both 
homogeneous and inhomogeneous deformations are pos- 
sible. In a general case, the elasticity-driven interaction 
between vortices include contributions due to both homo- 
geneous and inhomogeneous deformations. In the elasti- 
cally isotropic case the contribution due to the inhomoge- 
neous deformations vanishes. But in any case, the order 
of magnitude of the total interaction coincides with that 
of the contribution which is due to the homogeneous de- 
formations. 

The consideration of homogeneous deformations pro- 
vides us, in addition, an useful technical trick. It is based 
on to evaluate first the VL energy in the case in which 
the superconductor is elastically isotropic and its shear 
modulus (i is infinite. In this case, the calculations are 
free of approximations and are almost trivial. If fi = oo, 
the only elastic degree of freedom of the sample is its ho- 
mogeneous dilatation. Therefore, using already known 
formulas for the VL energy and taking into account the 
dependence on the dilatation of the corresponding coef- 
ficients, the elastic contribution can be easily computed. 
As we shall see, any isotropic case can be reproduced from 
this [i = oo one. Moreover, the previously reported re- 
sults can be easily checked by evaluating them for fj, = oo 
and comparing them with those obtained considering this 
case from the beginning. 

Let us mention also that we use the Fourier method 
when computing the VL energy in the elastically 
anisotropic case. This method permits to satisfy quite 
easily the boundary conditions of the elastic problem 
which correspond with those that take place in real ex- 
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periments. Thus, one avoids to reproduce spurious effects 
that a lack of attention to these conditions might give. 
One of such effect is, for instance, the sample form depen- 
dence of the elasticity-driven interaction between vortices 
(the same group of authors reported this dependence in 
Refill but not in Ref-12). 



II. ON THE ELASTIC EFFECTS WITHIN THE 
LONDON LIMIT 

When studying the influence of the elasticity on the 
vortex properties, many authors use an assumption which 
might seem quite natural (see, e.g., Refsi2*£*i2iii*l£). It 
consists of using the "London approximation" introduced 
by Abrikosov in Refill (see also Refill). However, the 
essence of this approximation could easily be misinter- 
preted. As it is frequently commented, within the Lon- 
don approximation the order parameter modulus varies 
significantly inside of the vortex cores only. Since the 
spontaneous deformation associated with the supercon- 
ductivity is proportional to square of the order parame- 
ter modulus, it seemed natural that only the core regions 
0° iS are essential sources of stresses. It is just what 
is assumed in Refsi2*£ii2*ii*i£. However, one has to bear 
in mind that supercurrents also produce an elastic effect 
because they diminish the value of the order parameter 
modulus. Locally this diminishing is small. But since the 
supercurrents occupy a very broad region (p < Al), their 
effect might be comparable and even more important, as 
virtually proves to be, than that of the cores. 

To make this point more clear, let us recall how the 
vortex self-energy per unit length Eq is calculated within 
the London limit Within this limit one assumes that, 
when calculating the supervelocity v s from the Ginzburg- 
Landau equations, the density of superconducting elec- 
trons (the square of the order parameter modulus f 2 ) 
is constant in the corresponding equation. This makes 
it possible to find out explicitly the spatial distribution 
of the supervelocity. After doing so, one can follow two 
different ways: 

(a) Following de Gennes?i& the vortex self-energy is 
presented as a sum of the magnetic field energy and 
the kinetic energy of the superconducting electrons: 



calculated from the exact formula 



e = J {H 2 + f 2 v 2 ) d 2 p 



(1) 



[we use here the reduced units, see Refil^l^, which 
are analogous of those defined in Eqs. (fTT)f> (see be- 
low)] . Integration is carried out taking into account 
the already found supervelocity, and considering 
that the density of superconducting electrons is 
constant. This approximation is justified by virtue 
of the high value of k: f 2 diminishes significantly 
only at p < £, whereas v 2 does at p > 

(b) Following Abrikosov (1^1^ the vortex self energy is 



£o = 



IF 



(1 - f 



d 2 p 



(2) 



(as before, we use here dimensionless quantities). 
The principal part of this integral arises from the 
second term, and it is associated with distances 
much larger than £. In other words, those varia- 
tions of / that takes place out of the vortex core 
are now essential. 

As we see, to assume that within the London approxima- 
tion / is constant out of the vortex cores is not always 
correct. But, as we have pointed out, this is just the 
assumption that unfortunately many authors made. For 
example, when studying the interaction between vortices 
and lattice defects, Miyahara et al~ considered integrals 
which arc similar to Eq. J3J) but, at the same time, ne- 
glected all the spatial variations of / at p > £. 

It is quite surprising that this assumption has not been 
critically revised up to now, especially by noting that, in 
principle, the importance of the out-of-core region for 
the elastic effects could be understood since long ago. 
Galaikc 4 considered the interaction between vortices and 
dislocation-induced strains. He found that this interac- 
tion depends not only on £, but also on A^. However, he 
did not comment Reff2> and discussed neither the vortex- 
induced strain nor the strain-induced interaction between 
vortices. ReffS* is a recent example in which the out-of- 
core region is taken into account when studying an elas- 
ticity related problem: the structure of a superconduct- 
ing vortex pinned by a screw dislocation. 



III. ONE SINGLE VORTEX 
A. Vortex-induced strain 

Let us proceed with the calculation of the strain field 
induced by one single vortex. When doing so, we shall 
account for all the spatial variations, core and non-core 
ones, that are associated with the vortex. 

The free energy can be presented as 



F = F l +F 2 = ^ J T 2 )dv, 



where v is the volume of the system, and 



-ihV-— A)* 



(3) 

(4a) 
(4b) 



Here and below, summation over double indices is im- 
plied. 
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The equations of equilibrium rea< 



115.17 



omit primes in the new quantities) 



a + b\^\ 2 + aqua + — ( - iKV— —A 



V x H 



1 / .„ 2e < - 

Am \ c J 

2i ' c 

Xijklim) + a ij (\^\ 2 ) = 0, 



* = 0, (5a) 



*| 2 A 







0. 



, (5b) 
(5c) 
(5d) 



where (...) means volume average. We shall look for the 
solution of these equations for the case of a single vortex. 
The z-axis of the coordinate frame we choose is parallel 
to the vortex. The crystal frame is obtained from this 
coordinate frame by rotation. 

It is clear that far enough from the vortex both the 
order parameter and the strain tensor tend to constant 
values; say $> s and ufj respectively. Assuming that 



(1*1 



\*S? S \ i the equations of equilibrium reduce to 



+ 6|* 8 | 2 + Q!ijUjj = °> 

x ijk iu s M + I 2 = o. 



In consequence: 



-a/b* 



u tj = atXhiXiju/b*, 



(6a) 
(6b) 



(7) 

(8) 



where b* = b-a l3 a k i\^i {\> u is given by \> u \ijm> = 
Skk'Sw). These values are just what one obtains in the 
superconducting phase. 

Putting Uij = + u^j , we can rewrite the equation of 
equilibrium l|5a|l as 



1 - 



1*1 



a ij u ij 
Ws\ 2 b 



* = 0, (9) 



where £ 2 = h 2 /(4m|* s \ 2 b). It is convenient to introduce 
the following notation: 



mc 



87re|* 



T7j 



H' 



H 

V2H C 



a 



A' = 



A' 



he 

2\/2e£A L 
_ _r_ 

A 

V2H C \ L ' 
A 



(10) 



Thus, the equations of equilibrium can be written as (we 



{l-vl-a ijU l)f-f = -K-^f, 

V x H = v s / 2 , 

K]kl(ukl) + tty(/ 2 ) = 0, 

d 

0^7 (Xtjkiu k i+a tj f 2 ) = 0, 



(11a) 
(lib) 
(11c) 

(lid) 



where the order parameter has been expressed as ^ = 
/e IX , with v s = k _1 Vx — A the above mentioned super- 
velocity. Here k — A^ /£ represents the Ginzburg-Landau 
parameter in our case, which does not differ substantially 
from the conventional one (b* ~ b). 

The spatial distribution of the supervelocity v s can be 
obtained from Eq. Ijllbfl by assuming that / is constant 
there, i.e. within the London limit. Thus one finds that 
v s = K~ 1 Ki(p), where K\ is the MacDonald function 
(see, e.g., Ref^). 

In Eq. I|llal) , the term with results to be of order 
a 2 because of Eqs. I|llc|l and (|lld(l . Since a is small^ 
the vortex-induced strain can be calculated to the lowest 
order in a neglecting the changes in / that the term 
in Eq. (|llafl induces. In other words, f 2 in Eqs. (|llc|) 
and l|lld|) can be taken as the solution of Eq. (|llafl with 
a = 0. This solution can be written as f 2 = 1 — h, 
where h represents the vortex contribution. Using the 
same approximation that in Ref^ we have 



Hp) 



v 2 s (p) 

i - cm 2 , 



P > K 
P -C K~ 



(12) 



where C is a constant of order unity. 

We present the vortex induced strain asi^ 



ij = e ij + \ fe w j(q) + * u i(q)] e 



q#0 



iA f 

€i 3 + g^2 / ^ fe u J'(q) + 1jUi(q)} ■ 



e iqp , (13) 



where A represents the section of the sample in perpen- 
dicular to the vortex (recall that we have split the total 
strain into itf • + u?-). Here ey accounts for the homoge- 
neous deformations that the vortex induces, and Ui(q) is 
the i-th component of the displacement vector in Fourier 
space. Thus, Eqs. (|llcl) and (|lld() can be written as 



Kjkim - otij (h) = 0, (14a) 
G?tt(q)u fc (q) + iSi(q)h(q) = 0. (14b) 

where Sj(q) = a^qj, G^(q) = XijkiQjQh and h(o) is tnc 
Fourier transform of the function l|12|). For the strain 
field we have: 



= a ki\ 3 \i{h), 



i(q) = ~iS k {q)Gki(q)h(q). 



(15a) 
(15b) 
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When calculating the strain tensor at a fixed distance p 
from the vortex [see Eq. (|13fl ], the inhomogeneous defor- 
mations are mainly given by those terms with q <C p . 
So the main contribution at p 1 arises from g -C 1. 
For these small q's, the function /i(q) can be split into 
core and non-core contributions: 

h core (q)~±- P f 27I (p-K 2 p 3 )e-^ cose dpde 
A Jo Jo 

= aJ ^-^ Jo ^^2^' ( 16a ) 



''lion — core 



(q) 



i 



2n 



Ak~ j k -i Jo 
2tt Z" 1 J (gp) 



p - l e- iqpcose dpd6 



dp = 



2tt 



<?«1 Ak 2 



hiK (16b) 



(here we have used the asymptotic form of v s ~ I /(up) 
for < p < 1, see Ref^). 

As a result, at p 3> 1 the strain tensor can be written 

as 



'/ 



OtklK A h 



ijkl 



A 



rf 2 q 



gA(q)G fcj (q)e iq -'' 



ijkl 



A 



. / Qij(6q)d0q 
P JO 



(17) 



where V — J h{p)d 2 p — n(l + 41iik)/(2k 2 ), and 0jj is 
a tensor which depends only on the angle # q ( q • p = 
qp cos 9q). If the sample is large enough the first term 
in Eq. i|17|l can be neglected. But we retain it because, 
when dealing with the strain-induced interaction (see be- 
low), its contribution becomes significant (this fact is well 
known in the theory of point defects, see e.g. Refsi2£). 
Note that the non-core contribution to r], i.e. the log- 
arithmic term, could also be obtained from the well 
known expression of the vortex self-energy: according to 
Abrikosov^ e ~\ J (I - f A )d 2 p ~ 2tt h{p)pdp = 

27TK~ 2 In K. 

Kogan et aliik obtained a similar expression for vortex- 
induced strain considering an infinite medium. In such 
a case, the first term of Eq. l|17|) vanishes at all. But 
the main difference between Eq. I|17f) an d the expression 
reported by Kogan et al. 11 resides in the corresponding 
values of rj. Assuming that only the vortex core induces 
strain, Kogan et al. reported a value tt/k 2 . So they over- 
looked the logarithmic term in 77 = tt(1 + 41n/-t)/(2K 2 ) 
that arises because of non-core contributions. This im- 
plies that in the case of high-K superconductors, Kogan 
et al. strongly underestimated the vortex-induced strain. 



B. Elasticity-driven interaction between vortices: 
Qualitative estimations 

Let us now estimate the interaction energy of a VL 
which is associated with the strains that the vortices in- 



duce. As we have pointed out before, the inhomogeneous 
part of these strains have been previously reported but 
neglecting non-core contributions (see, e.g., Refiii). If 
the distance between vortices is much longer than Al, to 
take into account these non-core contributions reduces to 
modify the strains by a factor. In consequence, the in- 
teraction energy one obtains by taking into account both 
core and non-core contributions coincides, up to the cor- 
responding factor, with previously reported ones. Kogan 
et a/.^i for instance, evaluated the interaction energy of 
a VL by summing up all pairwise contributions. Mod- 
ifying this interaction energy by including the non-core 
contributions, one can see that 



F 



(nh) (1 + 41n^) 2 AK 
k 2 K 



B 



(18) 



Here AK/K stands for the order of magnitude of the 
relative change in the elastic moduli due to the normal- 
superconducting transition, and B represents the mag- 
netic induction. 

The interaction between vortices that arise due to the 
homogeneous strains can be easily estimated as follows. 
It is clear that N vortices will induce a total (homoge- 
neous) strain Ne, where e is given by Eq. (|15af) . if the 
distance between them is large enough. When substitut- 
ing this strain in the corresponding terms of VL energy: 
-a(Ne)(N(h)) + \{Ne) 2 /2, one obtains -n 2 r] 2 a 2 /(2\), 
where n = N/A is the vortex density. This is precisely 
the interaction term that we are looking for. Taking into 
account that the vortex density is n = kB/(2it), and 
a 2 /\ = AK/K [recall that we are using the dimension- 
less units defined in Eq. IjlOl) ]: this interaction can be 
estimated as 



F 



(h) (1 + 41n«) 2 AK 
k 2 K 



B . 



(19) 



As we see, the order of magnitude of both interaction 
terms Eqs. (I18fl and (|19fl coincide. Consequently, either 
of them give us an estimate of the order of magnitude of 
the total interaction energy. 



IV. VORTEX LATTICE: ELASTIC ALLY 
ISOTROPIC MEDIUM 



It is convenient to begin the treatment of VL's consid- 
ering the case elastically isotropic superconductors. In 
this case, the elastic contribution to the VL energy can 
be obtained, without any new approximation, from al- 
ready known formulas for this VL energy. Such formu- 
las are available for the regions H w H c i, H w H C 2 
(Refsii^i^) and for intermediate fields H c i < H < H c2 
(RefsiiSiSi). They reasonably match at the boundaries 
of the corresponding regions. This permits us to study 
the elastic effects in isotropic superconductors with the 
same accuracy. We begin with the case p, = 00 where the 
calculations are elemental. 
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A. Infinite shear modulus 

The only elastic degree of freedom of a system which 
shear modulus is infinite is its homogeneous dilatation. If 
the system is not clamped this homogeneous dilatation, 
say u, must be understood as a variational parameter. In 
the free energy @, this variational parameter modifies 
the coefficient of the term Ivf^ 2 , which can be rewritten 
as a(u) = a + au. 

Let us fix the parameter u for a while, i.e., let us con- 
sider momentaneously a clamped sample. Thus, after 
minimizing with respect to all degrees of freedom except- 
ing u, the free energy of the VL with respect to that of 
the superconducting state can be written as a sum of two 
terms: a u-dependent VL energy via the coefficient a(u), 
and the elastic energy. It is (see, e.g., Reff 2 ^ and the 
references therein) 



F = F VL (u) + ^-u 2 



(20) 



where the Fvl has the form 
' BHd 



47T 



i 



1 



B + BH C \ — j- 



B 2 



1 



Cgc2- 

+ (2k 2 



B? 



(I), 



(I-II), 



(II). 



(21) 



in the corresponding regions of magnetic fields defined 
as (I): H ra Hoi, (I— II): H eX < H < H c2 , and (II): 
H « ff c2 . Here /3 A = (* 4 )/(* 2 ) 2 = 1.16 for a trian- 
gular VL, and 2 hi!/ = 2( 7 - 1) + ln[V3/(87r)], where 
7(= 0.57772 . . . ) is the Euler's constant. The magnetic 
induction B and the distance between vortices d are such 
that B = 2<j)Q I (s/3d 2 ) in a triangular VL, where </>o is the 
flux quantum. 

The magnetic induction as a function of the magnetic 
field is given by^ 



B 



20o 



In 



H — H c i 



30o 



Ait\ 2 L (H-H ol ) 



0o f 



In 



iTT\j{H-H cl ) 



00 



+ 7 



(I), 

} (I-II) 



tt _ H c2 -H 



(II), 



(22) 



where 7 = 2(1 — 7). Let us recall that in high-K super- 
conductors one has the following relationships (see, e.g. 
Refsi^): 



Zk 



4tt H c2 
V3 B 



8tt k 2 H c1 
y/Shxn B 



(23) 
(24) 



One can check that at the boundaries of the magnetic 
field regions, the expressions in Eq. 1)21(1 match one each 
other with a reasonably accuracy: 

o For H w (iH cl (Ci > 1) one has B w dH cl ac- 
cording to the expressions (I-II) and (II) in Eq. 
J32J). Therefore, taking into account the relation- 
ship J23J, one can obtain the free energy as 



F 



VL 



An 



SIT I « 



In 



8it(vk) 2 /(V3(i ln/c) 



m k 



(I), 

(I-II) 

(25) 



o For 7J = ( 2 H c2 (C2 < 1) one nas ^ — QiBci accord- 
ing to the expressions (I-II) and (II) in Eq. i|22|) 
and the relationship J23J. Therefore, taking into 
account the relationship l|24l) . one can obtain the 
free energy as 



& H c2 I ^ + 



Ft 



VL 



H 2 2 
~5F 



ln[47n/ 2 /(V3C 2 )] 



2k z 



C 2 - 



d-c 2 ) 2 

2k 2 



(I-II), 



(II). 



(26) 



The critical magnetic fields entering all above expres- 
sions are w-dependent magnitudes: 



Hci(u) = — =- 
\J2n 

H c2 (u) = V2kH c (u) = H° 2 + H' c 



H c (u) = H° cl + H' cl u, (27) 

(28) 



where H c (u) = 2a{u)y^F/b = H° + H' c u (with H° = 
2ay/ir/b and H' c = 2a^J 7r/6). In consequence, the ratio 
d/£ is also a w-dependent magnitude [see Ea. lf^Ijl ]. from 
which one can write the coherence length £ as £ = £° + 
£'u. Mention that the Ginzburg-Landau parameter k is 
independent of u because it does not depend explicitly 
on the coefficient a(u). 

Let us now proceed to minimize the free energy <|20[) 
with respect to u, i.e., to take into account that the sam- 
ple is in fact undamped. After doing so, we obtain 



' ^{2BH° cl -5 l B 2 



F = < 8¥ 



1 

8¥ 



B 2 + BH o ci H^mi _ , s 



In k 



(I), 



(HI), 



B 2 — 



(H° c 



B) 



1 + (2k Z -1)Pa~ Pe 



(II), 



(29) 



6 



where 



Si-ii 



In 2 k AK 

2k 2 K ' 

[l + 21nH/£°)] 2 i^ 2 



ln 2 (d/£°) AK 



AK 
~K 



2V2H' C 2 B 



4k 2 



K 



(3 e = 2k 2 — 



(30) 
(31) 

(32) 



Here it has been taken into account that H' 2 /(AttK) = 
AK/K is the relative change in the bulk modulus due 
to the normal-superconducting transition. Because this 
relative change is usually AK/K <C 1, the expression for 
the region (II) in Eq. (|29[) can be written as 



F ~ — 

8tt 



B 2 



B) 



1 + (2k 2 - 1)0 A 



where 



[1 + (2k 2 - l)p A ] 



Sn(H c c 



1 AK 
~2^~K' 



B) 



(33) 



(34) 



In all above expressions for the free energy, one can 
identify a term 



SB 2 



(35) 



which describes an attractive interaction between vor- 
tices. The coefficient 5 is given by Eq. (J3DJ, (J31J or 
depending on the magnetic field region one considers. It 
can be presented as 



[( + \n(d/C)} 2 AK 



2k 2 



K 



(36) 



taking into account that the ratio must be replaced 
by k if d > Al, and £ is a constant of order of unity (see 
Fig. P). 

According to what we have seen in the preceding sec- 
tion, the logarithmic contribution to the coefficient S is 
due to non-core effects. These effects have been over- 
looked until now. As we show in FigQ] the neglec- 
tion of these non-core effects leads to underestimate the 
elasticity-driven interaction between vortices. And by 
virtue of the high value of k, such a underestimation is 
quite significant in almost all the mixed state. 



100 r 



TTTT] 1 TT-I 

r 
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1 I 1 1 1 1 1 1 

(II) 1 
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FIG. 1: Log-log plot of the coefficient 5 of the attraction term 
~ -SB 2 of the free energy as a function of the magnetic induc- 
tion, taking into account (solid line) and neglecting (dashed 
line) non-core contributions. The regions indicated as (I), (I- 
II) and (II) (see text) are the corresponding ones for k ~ 100 
(note that H c \ ~ 10 _4 H C 2 in this case). 



where K4/3 = K + 4/i/3. The first term of this expres- 
sion renormalizes the coefficient b of Eq. 0]) . This renor- 
malization disappears in the limit fi — * 00. The second 
term makes that the free energy becomes into a non-local 
functional. This non-locality remains as long as the shear 
modulus does not vanish. 

Working out from this functional one could recover 
the free energy l|29l) by ascribing the coefficients in Eq. 
(|37|l the values they assume in the case /i = 00, i.e. 
a 2 /(2JQ/ 3 ) = and 4^/(3X4/3) = 1. This consideration 
does not change the functional form of Eq. (|37|l (only co- 
efficients change). So one concludes that the free energy 
density of any isotropic type-II superconductor has the 
form of Eq. (|29|l with the corresponding renormalized 
constants: 



(a 2 /K) 



b-a 2 /K 4/3 (38) 
(a 2 /^)[4/i/(3^ 4/3 )]. (39) 



Note that, because the resulting coefficient 5 in Eq. i|35|) 
vanishes if /i = 0, it can be said that the elasticity-driven 
interaction between vortices is associated with the solid- 
state elasticity. 



B. Finite shear moduli 

It is quite straightforward to extend the results that we 
have obtained for the (J, = 00 case, to the most general 
isotropic one. Note that minimizing the free energy @ 
with respect to all elastic degrees of freedom one obtains 

^ = -4l <l * lV ^4l <l *' 2>2 ' (37) 



C. Comparison with previously reported results 

Let us start this section by comparing our results with 
those reported in Refsii2*ii. In these references, interme- 
diate fields far from T c are considered. Mention that al- 
though strictly speaking the Ginzburg-Landau approach 
that we use is not valid far from T c , it still gives correctly 
the orders of magnitude. So the comparison still makes 
sense. Mention also that in Refs>i2iii the homogeneous 
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part of the strains are omitted. In RefiiS this omission is 
mentioned explicitly, while in Ref^i it follows from the 
fact that they consider infinite samples when calculating 
the interaction between vortex pairs. Therefore, as we 
argued in Sec. IIII Bl we can only compare the order of 
magnitude (see below for a more detailed comparison). 
Such a comparison reveals that, as a result of the neglec- 
tion of the non-core contributions to the interaction en- 
ergy, this energy is notably underestimated in Refs>i2iii 
through the most part of the mixed state (see Fig. 0. 
Such a underestimation is at least by factor ~ In 2 k close 
to H cl . 

In Refpi^, treating the case H H c2 , both homogenous 
and inhomogeneous strains are seemingly taken into ac- 
count. In accordance with Eq. (28) of this reference, the 
free energy in the isotropic case should be of the form 



F = — 



where 



B 



1 + (2k 2 - l)0 A - 4k 2 /3 2 
[1 + (2k 2 - 1)0 A + 4k 2 /? 2 ] 



;(H° c2 -B) 



0A 



(40) 



(41) 



Because of the smallness of 02, this expression can be 
approximated to 



8tt 



B 2 - 



(H° 2 - B) 



1 



(2k 2 



-5 2 {H° 2 -Bf 



(42) 



where 8 2 = -I2n 2 2 /{1 + (2k 2 - 1)0 A } 2 . 

At first glance, one could think that the last term of 
this expression and the last one of Eq. IM/il) differ only 
in a numerical factor. So the elasticity-driven interaction 
between vortices are reasonably well reproduced by ei- 
ther of them. However, a deeper inspection of Eq. i|42|) 
reveals that it is erroneous: the coefficient S 2 (i) van- 
ishes if fi — oo and (ii) remains finite if /i = 0. Being 
the strain-driven interaction between vortices due to the 
specific features of the solid-state elasticity, the results 
one obtains from Ref^£ in the above mentioned limiting 
cases cannot be correct. This motivates us to reconsider 
the problem treated in Refii 2 - (Sec. IV Bl below 1 ) . 



Because $ satisfies the equation (|5all , this expression can 
be written as 

1 r / 9p \ 2 

— /W -ihV A) ^dv 

4m J \ c 1 

(a|*| 2 + 6|*| 4 + a ijUij \V\ 2 ) dv. (44) 

As a result, the free energy can be presented as 
1 ffH 2 1 



F = 



I \ N/r ~ 2"'*' 4 + ^vkiUijUH ) dv. (45) 



This expression generalizes the Abrikosov' one [see Eq. 

by taking into account the elastic degrees of freedom. 
In Refi^ it has been reported a similar expression that, 
however, is erroneous (see Sec. IV Bl below 1 ) . 



A. H < H c2 

When treating the fields far from H c2 , it is convenient 
to express I*! 2 = |* s | 2 - h, where |\I> S | 2 = -a/b* [see 
Eq. JJJl] and h now represents the VL contribution. In 
addition, we express the strain tensor as ity = uf^ + u"-, 



where < = -a k i\L |*J 2 and 



~ [qiS k (q)G kj (q) + qj S k (q)G H (q)] h(q)e lc 



q#0 



(46) 



Thus, the equations of equilibrium (|14a|) and (|14bjl are 
satisfied. Substituting these expressions for l^l 2 and Uij 
into Eq. Q45|). one finds that 



i^6'(q)|Mq)| 2 , 



(47) 



where F s = -6*|* s | 4 /2 and 



V. VORTEX LATTICE: ELASTIC ALLY 
ANISOTROPIC MEDIUM 



, U^A-L (q=0), 

W U(q)^(q)G^(q) (q^O). 1 ' 



Let us reconsider the free energy Q and, as usual (see 
Refii^), integrate by parts term with V^*. Thus, after 
using the Gauss' theorem and the boundary condition 
n • ( — ihV — ^rA)^| s = (n is the unit vector of the 
normal to the surface £), one finds that: 



1 

4m 



2e , 
ihV A j M> 



dv 



1 

4m 



2e \ 2 
— A) Vdv. 

c 



(43) 



Note that, for q ^ 0, the function b'(q) only depends on 
the angle. 

When calculating h, we can retain the lowest order 
terms in a, i.e. we can take h ~ ho + hi where ho is 
given by solving Eq. (5aJ) with a = 0, and hi represents 
the correction to this solution due to the term onjV^A in 



Eq. 1(53 



a ij u ij ( ou t °f cores), 
(inside cores). 



(49) 
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The absence of correction inside the cores follows from 
the fact that, in these regions, Eq. I(5a() can be linearized 
because of the smallness of the order parameter modulus 
(see, e.g., Ref^). Let us remark, even doing so, the 
vortex cores are taken into account: they act as strain 
sources. 

Substituting these expressions in Eq. (|47|l , and retain- 
ing the lowest order terms, we obtain 



F =F s + b*\H> s \ 2 ((h Q ) + 



(H 2 ) 
8tt 



;(h 2 ) 



& £ V-q)^i (q) + 1 E 5 ' 1 2 ' 



(50) 



where £ = 1 + b' Q /b. 

The most important contribution to the two last terms 
in Eq. f5T)| arises from q < At these q's, the 

function /ii(q) can be calculated by taking bh± ~ a^u^ 
in all the regions [see Eq. 149(1 ] . so 



-&£>o(-q)Ai(q) ~ -J2b'(q)\h (q) 



(51) 



As a result, the free energy l(57)|) is 
F~F s + b*\* s \ 2 ((h ) + 



(H 2 



8tt 



(52) 



The last term of this expression represents the strain- 
induced contribution to the VL energy. The term with 
q = is associated with the homogeneous strains. The 
elastic constants enter this term through an invariant 
combination [see Eq. 1(48(1 ] . so it does not depend on the 
orientation of the VL with respect to the crystal axes. 
This dependence arises from the terms with q ^ 0. 

Let us mention that this formula demonstrates that the 
elasticity-driven interaction between vortices does not de- 
pend on the sample form, unlike to the statement made 
in Ref^-. Indeed, such a dependence would mean that 
contribution to the sum from the region of small q's is es- 
sential and comparable with the contribution of the rest 
of the sum. But the function /i (p) — (h ) is a periodic 
function defined in a finite volume (neglecting the near- 
of-the surface distortions). Its Fourier spectrum does not 
contain small g's but has maxima at the non-zero recip- 
rocal lattice vectors. The form and the size of the sample 
is reflected in the form and the width of these maxima 
and nowhere else. In fact, the sums over q's can be re- 
placed by sums over the reciprocal lattice vectors of the 
VL Q. Putting 



ho(p) = ^2h {p- Pi) 



(53) 



where p i represent the vortex positions, one finds that 
h ^ = 1^ / ~ h ^P - Pi^'^P = nh (q), (54) 



where n is the vortex density, A represents the section of 
the sample in perpendicular to the VL, and 



fro(q) 




(q=Q), 
(otherwise), 



(55) 



with Q any of the reciprocal lattice vectors (note that 



Eq. 1(52(1 can be written as 



nSqQ). As a result, the last term in 



F r] = - — 



b'(0)h 2 (0) 



£&'(Q)MQ)I S 

Q#0 



(56) 



Let us emphasize that with this expression, one takes 
into account that both core and non-core regions act as 
strain sources. It can be straightforwardly illustrated 
close to H cl . Here, due to the large separation between 
vortices, the function h in Eq. ((53(1 practically coincides 
with that associated with one single vortex [see Eq. (|12|l ] . 
In consequence the function /io(Q) varies slowly up to 
Q « £ _1 and then rapidly drops to zero. So in Eq. 1(561) 
it can be taken as /io(Q) — ^o(O), which naturally split 
into core and non-core contributions [see Eqs. JTBJ]: 



ho(0) =^(l+41n K )|1< s | 2 e72, 



limiting the sum over Q's up to <5 max ^ S,^ 1 - 

Let us calculate F c \ explicitly in the isotropic case. In 
this case one has ctij = aSij and Xijki = (K — jjj)8ij8ki + 
^(Sik^ji + SuSjk), where K and /i are the bulk and the 
shear modulus respectively. In consequence, 



(57) 



& '(Q) = ^ 2 



and Eq. l(5T)|> yields 



F, 




4/x/3) 



(Q = o), 
(Q^o), 



(58) 




a 2 hl(0) Mfin 2 r 



(59) 



where the sum over discrete Q's has been replaced by in- 
tegration Q^q ~ n^ 1 j gPQ). The term oc n represents a 
renormalization of the vortex self-energy, while the term 
cx n 2 is the elasticity-driven interaction. 

In the anisotropic case the sum over Q's in Eq. I(56|) 
also yields a term oc n, which renormalizes the vortex 
self-energy, and a term oc n 2 which contributes to the 
elasticity-driven interaction between vortices. In this 
case, a dependence on the orientation of the VL with 
respect to the crystal axes is implicit in these two terms. 
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Taking into account that |* s | 4 6' ~ H 2 (AK/K) and 
nH c ^ 2 ~ B/k, the elasticity-driven interaction can be 
estimated as 



(l + 41n«) 2 AK p2 

Fint ^ ic B ■ 



(60) 



As we see, its order of magnitude coincides with that we 
obtained in Sec. IIII Bl from qualitative estimations, as 
well as with that of the exact results that we obtained in 
Sec. II VI for the isotropic case. Let us mention that omit- 
ting In k in these expressions, i.e. omitting the non-core 
contribution to the vortex-induced strain, they reproduce 
the previously reported results^i 



B. H « H c 



Here Hq is a constant and u> is the squared modulus of the 
function which is solution of the linearized equations 
(I5ajl and l|5b|) (\P = ^/we ix ). Bearing in mind that the 
magnetic induction is B — (H z ) = Hq — (lo)/(2k), from 
Eqs. l|6*T|) one can also obtain the following relationship: 



2k(k-B) 

0A~f3e 



(62) 



(2k 2 



where f3 e = —2n aij{uijUi)/{u)) and Pa = 
1)0 A , with p A = (cc 2 )/» 2 . 

Bearing in mind that from the Abrikosov identity l|61a|) 



it follows that (H 2 ) = B 2 + ((uj 2 ^ 
energy (|45|) can be rewritten as 



(w) 2 )/(4k 2 ), the free 



Let us now consider the VL's near H c i- When doing 
so, it is convenient to use the conventional dimension- 
less units instead of those defined in Eqs. (|l(Jfl . These 
conventional units can be obtained from Eqs. (|10J) by 
replacing \^ s \ 2 with —a/b. 

Following Kogan24 one can easily obtain, now from 
the equations of equilibrium © , the so-called Abrikosov 
identities (see also Refj^i^S^) . In presence of strain they 
read 



k - H 



H z = Hq 



2k 



(61a) 



1 - 2k' 
2k 2 



(uj 2 ) - on {uijuj) = 0. (61b) 



F = B 



1 



4k2 M + 2 X vki(u t3 u k i). (63) 



When calculating the strain tensor u%$ (now the total 
strain), we must take into account that l^l 2 = u>. Thus, 
expressing Uij in the form (|13|l . one finds that 



ijkl 



Ui(q) = iS j (q)G ji (q)u)(q), 



(64a) 



where w(q) represents the Fourier transform of the func- 
tion uj. In consequence, the last term of Eq. (|63|l is 



Kjkl{UijUkl) 



1 r 

2 



Xijkieijeki + G i/(qK(q)%(-q) 

Ui(q)u>(-q) 

a ij e lj (uj) + i^Sj(q)u i (q)a;(-q) = -~aij{uiju) = -^(oj) 2 . 

q^O 



(65) 



r 



As a result, the free energy Ij63(l is 



Ak 2 



0A~Pe 



(66) 



Let us mention that the form of this expression for 
the free energy differs substantially from that reported 
by Miranovic et al. in Ref. 12 . The elasticity-driven in- 
teraction term that one obtains from Eq. I|66|l is several 
times smaller than the corresponding one in Refii^. The 
validity of Eq. I|66|l can be checked by noting that it 
reproduces the isotropic case [see expression (II) in Eq. 
(|29p]. In contrast, the expression reported in Refii^ does 



not (see Sec. IIV C|) . The reason is that is it obtained 
from an expression analogous to Eq. (|45|l . but erroneous 
[Eq. (20) of Ref. 12 ]. It reads 



F 



H 2 b 1 \ 

o ol*| 4 + ctijUij |*| 2 + -A ,,v; a,, in,, I dv, 

07T Z Z I 

(67) 



One can see here that the term aijUij\^\ 2 is taken into 
account twice: one explicitly and another implicitly in 
the term — £>|<I'| 2 /2 which arises as result of the integra- 
tion by parts showed at the beginning of this Section. 
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VI. CONCLUDING REMARKS 

We have revised the contribution to the VL energy 
which is due to the vortex- induced strains, showing 
that essential corrections in the previous calculations are 
needed. The most important one is connected with the 
fact that, in high-K superconductors, not only the vor- 
tex cores induce strains in a significant way. There also 
exists a significant contribution associated with the non- 
core regions which, in fact, the most important ones for 
the VL energies at low fields (H -C H C 2). As a result 
of the proper inclusion of all strain sources, the strength 
of the elasticity-driven interaction between vortices in- 
creases by a factor up to ~ In 2 k compared with the pre- 
viously reported ones. 

It is known since long ago that the observed correla- 
tions between VL's and crystal lattices in dirty supercon- 
ductors cannot be explained without the elasticity-driven 
interaction between vortices^S This interaction has been 
proved to be important in clean superconductors also. 
For example, the VL's observed in NbSe2 do not corre- 
spond to the minimum of the London energy. In Ref. 11 
Kogan et al. showed that, however, the difference in the 
London energies of the two possible competing structures 
is smaller than the difference in the energies of the corre- 
sponding elasticity-driven interactions. As we have men- 
tioned, Kogan et al. underestimated the elasticity-driven 
interaction between vortices because they assumed that 
only the vortex cores induce strain but, even doing so, 
they pointed out the importance of this interaction in 



NbSe2- This importance is increased as a result of the 
present work, what should be taken into account espe- 
cially in those cases in which previous estimates of the 
above mentioned differences concluded that it was that 
of the London energy the most important one. 

V 3 Si might provide an example in which the latter case 
takes place. In Ref. 2 ^« it was claimed that in VsSi the con- 
tribution to the VL energy which is due to the (under- 
estimated) elasticity-driven interactions between vortices 
can be neglected compared to that contribution due to 
the nonlocal corrections to the London energy. But bear 
in mind that (i) the order of magnitude of these two 
contributions is the same, as it was shown in Ref^i con- 
sidering the vortex cores as the only sources of strains, 
and (ii) the strength of the elasticity-driven interaction 
is considerably stronger than it was reported, as we have 
shown in this paper. So it is quite probable that in VaSi, 
as well as in other superconductors with large n, this 
elasticity-driven interaction between vortices is not only 
comparable, but even more important than the nonlocal 
corrections to the London energy. 

Acknowledgments 

We acknowledge L.S. Froufe for useful discussions. 
A.P.L. was supported from the ESF programme Vor- 
tex Matter in Superconductors at Extreme Scales and 
Conditions (VORTEX) and by Comunidad de Madrid 
(07N/005/2002). S.A.M. was supported from the Rus- 
sian Fund for Fundamental Research (Grant No. 00-02- 
17746). 



1 RW. Anderson, Phys. Rev. Lett. 9, 309 (1962). 

2 E.G. Kramer and C.L. Bauer, Philos. Mag. 15, 1189 
(1967). 

3 R. Labusch, Phys. Rev. 170, 470 (1968). 

4 V.P. Galaiko, Pris'ma Zh. Eksp. Teor. Fiz. 7, 294 (1968) 
[JETP Lett. 7, 230 (1968)]. 

5 K. Miyahara et. al., J. Phys. Soc. Jpn. 27, 290 (1969). 

6 A.M. Campbell and J.E. Evetts, Adv. Phys. 21, 199 
(1972), and the references therein. 

7 H. Kronmller and R. Schmucker, Phys. Status Solidi B 57, 
667 (1973); 74, 261 (1976). 

8 E.M. Chudnovsky, Phys. Rev. B 64, 212503 (2001). 

9 B. Obst, Phys. Status Solidi B 45, 467 (1971); J. Schel- 
ten, H. Ullmaier, and W. Schmatz, ibid 48, 649 (1971); U. 
Essmann, Physica 55, 83 (1971); H.F. Hess, CA. Murray 
and J.V. Waszzak, Phys. Rev. Lett. 69, 2138 (1992), Phys. 
Rev. B 50, 16528 (1994); CA. Bolle et al., Phys. Rev. Lett. 
71, 4039 (1993); P.L. Gammel et al., Phys. Rev. Lett. 72, 
278 (1994). 

10 H. Ullmaier, R. Zeller and P.H. Dederichs, Phys. Lett. 
44A, 331 (1973). 

11 VG. Kogan et al., Phys. Rev. B 51, 15 344 (1995). 

12 P. Miranovic, L. Dobrosavljevic-Grujic and VG. Kogan, 
Phys. Rev. B 52, 12 852 (1995). 

13 E. Simanek, Phys. Lett. A 154, 309 (1991); M. Coffey, 
Phys. Rev. B 49, 9774 (1994); J. Low Temp. Phys. 96, 81 



(1994); 97, 181 (1994). 
14 AA. Abrikosov, Zh. Eksp. Teor. Fiz. 32, 1442 (1957) [Sov. 
Phys.- JETP 5, 1174 (1957)]. 

A. A. Abrikosov, Introduction to the Theory of Metals 
(Noth-Holland, Amsterdam 1986). 

16 PG. de Gennes and J. Matricon, Rev. Mod. Phys. 36, 45 
(1964). 

17 L.D. Landau and E.M. Lifshitz, Theory of Elasticity (Perg- 
amon, New York, 1986). 

18 The corresponding parameter of smallness is AK/K ~ 
aijOtki\~j kl , where AK is the change in the bulk modu- 
lus K at the normal-superconducting transition. 

19 A.I. Larkin and SA. Pikin, Zh. Eksp. Teor. Fiz. 56, 1664 
(1969) [Sov. Phys. JETP 29, 891 (1969)]. 

20 AG. Khachaturyan, Theory of Structural Transformation 
in Solids (Wiley, New York, 1983); C. Teodosiu, Elastic 
Models of Crystal Defects ( Springer- Verlag, 1982). 

21 J. Friedel, PG. de Gennes and J. Matricon, Appl. Phys. 
Letters 2, 119 (1963); A.L. Fetter, Phys. Rev. 147, 153 
(1966). 

22 A.L. Fetter and P. Hohenberg, in Superconductivity, edited 
by R.D. Parks (Marcel Dekker, Inc., New York, 1969), Vol. 
2, pp. 817-923. 

23 Note that, in principle, d~ , A^ 1 , and £ _1 are characteristic 
q's of the functions /io,i(q) (d is the distance between vor- 
tices). Close to H c \ they are such that gT 1 <C A^ 1 <C 



11 



For intermediate fields, however, Xl is not a characteristic 25 V.G. Kogan and J.R. Clem, Phys. Rev. B 24, 2497 (1981). 

length of the VL because of the overlap between vortices. 26 V.G. Kogan et al, Phys. Rev. Lett 79, 741 (1997). 
So here we have dT 1 <C 
24 V.G. Kogan, J. Low Temp. Phys. 20, 103 (1975). 



